Parity-broken ground state for the spin-1 pyrochlore antiferromagnet 
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The ground-state properties of the spin-1 pyrochlore antiferromagnet are studied by applying the 
VBS-like tetrahedron-unit decomposition to the original spin system. The symmetrization required 
on every vertex is taken into account by introducing a ferromagnetic coupling. The pairwise effective 
Hamiltonian between the adjacent tetrahedrons is obtained by considering the next nearest neighbor 
and the third neighbor exchange interactions. We find that the transverse component of the spin 
chirality exhibits a long-range order, breaking the parity symmetry of the tetrahedral group, while 
the chirality itself is not broken. 

PACS numbers: 75.10.Jm, 75.40. Cx, 75.10.-b 
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The pyrochlore lattice, the network of the corner shar- 
ing tetrahedrons or a fcc-array of tetrahedrons, is a typ- 
ical example of a three-dimensional(3D) frustrated sys- 
tem, which is found in a number of materials, such as 
spinels, pyrochlores and C15-type Laves phase. To study 
effects of the geometric frustrations and the resultant 
enhanced spin fluctuations, the pyrochlore spin systems 
have been investigated intensively. For the spin-1/2 anti- 
ferromagnetic(AF) Heisenberg model, the quantum spin- 
liquid ground state is proposed based on the series expan- 
sion at finite temperature pj, but concerning the prop- 
erties of the low-lying spin-singlet states below the spin 
gap, a consistent picture has not emerged yet |g-|j- 

For the spin-1 case, numerical methods powerful for 
the spin-1/2 systems become difficult to obtain definitive 
results |ij and we need analytic approaches to solve the 
problem. We proposed the idea of the tetrahedron-unit 
decomposition of the pyrochlore lattice, which is a nat- 
ural generalization of the valence-bond-solid(VBS) state 
approach developed for the ID spin-1 systems M. In this 
approach the ground states of the fundamental unit, cor- 
responding to the valence bond for the VBS case, is the 
tetrahedron spin singlets which form the two-dimensional 
E representation of the tetrahedral group (Id). In conse- 
quence, the constructed VBS-type wavefunctions define 
the variational spin-singlets manifold with a macroscopic 
degeneracy, to which the ground state of the original 
problem may be continued adiabatically. In this situa- 
tion, it is essential to investigate how the ground state in 
the thermodynamic limit is stabilized from a collection 
of spin singlets by lifting the degeneracy. 

In the previous paper H, we examined the magneto- 
elastic couplings between the tetrahedron-singlets and 
the local lattice distortion of the E modes on the same 
tetrahedron as the source of lifting the degeneracy. As 
a result we have found that a Jahn- Teller mechanism 
driven by the non-magnetic spin degrees of freedom gives 
rise to the structural phase transition. The local lattice 
distortion is given by Q v mode of Tj. Actually the uni- 



form Q v distortion compressed along the c-axis is consis- 
tent with the cubic to tetragonal structural phase transi- 
tion without accompanying any magnetic order, observed 
in the spin-1 spinel-type antiferromagnets, Z11V2O4 and 
MgV2C>4. This scenario can explain, at least qualita- 
tively, the experimental results around the structural 
transition temperature. 

As another way of lifting the degeneracy, in this letter, 
we study the case where the inter-tetrahedron interac- 
tions, the next nearest neighbor and the third neighbor 
spin- spin interactions, are relevant. In particular, we will 
discuss the properties of the order parameters of the ob- 
tained ground state. Since the tetrahedron-singlets may 
be labeled by the spin chirality, a chiral ordered state is 
an interesting possibility. However, as long as physically 
reasonable interactions, such as spin exchange interaction 
or the dipole interaction, are considered, we find that it is 
not the chirality but its transverse component that shows 
the long-range order. 
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FIG. 1. The primitive cell of the pyrochlore lattice (a) and 
its projected view from the [lll]-axis (b). The small cubes in 
figure (a) are the guide for eyes with one-sixteenth volume of 
the unit cell. 



Following the ID VBS approach, we break up the orig- 



inal spin-1 objects (denoted by Si) into the two spin-1/2 
ones (si a and Sib) and introduce a ferromagnetic Hund 
couplings between them (—Jp), which serve as the sym- 
metrization. To project out the singlet sector of the 
composite spins completely it is necessary to take the 
Jp — > oo limit, however in nature a strong but finite Jf 
is sufficient. The pyrochlore lattice may be decomposed 
into the A and B types of bonds, which form the up- 
ward and downward tetrahedrons as shown in Fig. 0(b). 
Because of the symmetrization it is allowed to consider 
that the 1/2-spins with a (b) index construct the network 
of the A (B) bonds without loss of generality. In this 
decomposite-spin representation, the original AF Hciscn- 
berg model can be rewritten as follows in the strong fer- 
romagnetic coupling limit (Jf — > oo); 
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In the Jf — limit, the ground states of Eq. (|l|) are 
given by; 
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l*o) = Y[(a k \u) k + Pk\v)k) 



(2) 
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with arbitrary complex parameters au and (3k keeping 
\ a k\ 2 + \Pk\ 2 = 1] where TV is the number of sites and 
k specifies a tetrahedron. Here the orthonormal basis, 
{|u),|i>)} defined for a single tetrahedron, are the total 
spin-singlet states with s\ x + s 2x = s 3x + s ix = or 
1, respectively, (x = a,b) (see Fig. H). These tetrahe- 
dron singlets are the non-magnetic doublets belonging to 
the E representation of Td and the 2 JV / 2 -fold degenerate 
manifold defined by Eq.(||) is expected to be adiabati- 
cally continued to the low-energy sector of the original 
model. When representing the symmetry operation of 
T c i in this singlet subspace, {|w), \v)} real-basis diago- 
nalizes the parity operations with respect to the bonds 
vertical to the c-axis. On the other hand, the chirality 
basis M, which is defined by \R) — (\u) — i\v)) /V2 and 
\L) = \R)*, diagonalizes the four distinct C3 operations 
with the eigenvalues of oj = (—1 + voi) /2 and its com- 
plex conjugate uj*, respectively. 
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FIG. 2. The schematic representations of the tetrahe- 
dron-singlet states, |u) and \v). 



We proceed to the next step of lifting the 2 Ar / 2 -fold 
degeneracy of Eq. (g) . When considering the pure two- 
tetrahedron problem, we can easily see that the ferro- 
magnetic coupling Jp alone does not lift the degeneracy. 
Therefore the pairwise perturbation between the two ad- 
jacent tetrahedrons caused by Jf never fixes the original 
2 x 2-fold degeneracy. ^From the viewpoint of geometry, 
this means that the relative rotation of two tetrahedrons 
is not fixed by only Jp ■ The degeneracy is partly lifted by 
considering three tetrahedrons and eventually the degen- 
eracy will be lifted for the entire lattice. However energy 
scales of the lifting of degeneracy are expected to be small 
and hierarchical B . This problem is interesting but more 
academic since in nature there are other couplings to lift 
the degeneracy. In we have considered the local cou- 
pling with lattice distortion which leads to the structural 
phase transition. As another source of lifting the degen- 
eracy we can introduce longer range interactions which 
are relevant to lift the degeneracy already for a pair of 
tetrahedrons. 

In order to fix the relative rotation of the two adja- 
cent tetrahedrons, we include the next nearest neigh- 
bor (Ji) and the third neighbor interactions (J2). In 
spinels, J2 is expected to be important because of the 
existence of a superexchange-path through the oxygens 
on a single plane. The arrangement of these interactions 
are shown in Fig. |3|. By the second-order perturbation in 
— Jf, Jj, and J2 within the degenerate subspace spanned 
by Eq.(Q), we obtain the pairwise effective Hamiltonian 
between the tetrahedron-singlet states on the adjacent 
tetrahedrons as follows; 
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where the coupling is completely the same for the four 
distinct directions of pairs. In Eq.p), 



d = - — { (Jf - 4J 2 ) 2 + 2 (J F - 4^) 2 

+ 10(4J 2 -4J!) 2 }, (4) 



and the Pauli matrix ffc, describing the two-dimensional 
spin-singlet space on the kth tetrahedron, is defined by 
using the chirality basis as, 

t* = \R)(R\ - \L)(L\, 

t+ = \R)(L\, r- = \L)(R\. 

In fact, we can show that the universality class of the 
effective Hamiltonian is XX Z type by using a symmetry 
consideration as follows. First, the ±2-7r/3 rotation, C3 



and C~ , around the axis joining the centers of two tetra- 
hedrons leads to the conservation of their total chirality. 
Second, the \R) and \L) states are mutually conjugate by 
the time-reversal symmetry. Therefore we can conclude 
that the derived effective interaction has the XYZ-type 
symmetry of the pseudo r-spins and that the Ising term 
must vanish in the even-order perturbations. It should be 
noted that the second-order perturbation gives the lead- 
ing terms, since we start from the non-magnetic zeroth 
order states. 

In order to estimate the higher-order perturbation 
terms concerning Jp within the pairwise treatment, we 
consider the infinite Jp limit by considering the original 
spin-1 object on the sharing top vertex of the two tetra- 
hedrons. In this case, the effective Hamiltonian obtained 
by the first-order perturbation is given by, 
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( h - Ji) {rtj- b + t-jD + C 2 , (5) 



with C 2 = (J 2 + 2J0 /9, which is just 8/9 of the coeffi- 
cient of the first order term of Eq. ([|) in Jp. The r ma- 
trixes here are defined concerning the bottom triangles, 
which are related to the f matrix in Eq. ([|) through the 
relations like \R) = {\] R )\ [ v ) - \[ R )\ Tv» /V2, where 

ti?) is the doublet states with R chirality about the 
bottom triangle and | ]v) represents the spin state on 
the top vertex. It should be noted that such decomposi- 
tions cannot be extended consistently all over the lattice. 
Here we have shown it to illustrate the generic form of 
the pairwise effective interaction in the strong Jp limit. 

Comparison of Eq. (||) and Eq. (||) suggests that 
higher order effects of Jf only renormalize the strength 
of the pairwise interaction. Therefore, we assume that 
the effective Hamiltonian to describe the low energy part 
of the original Heisenberg model with next nearest neigh- 
bor (Ji) and third neighbor interactions (J 2 ) is given as 
follows, by using the spin-1/2 pseudospin operator, fu, 
defined for the fcth tetrahedron; 
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where J' — J 2 — J\ with a positive c. Since both of the 
upward (labeled by k a ) and downward (kb) tetrahedrons 
form the distinct fee lattice structure, Eq. (pf) is the 
spin-1/2 XY model on ZnS type bipartite lattice with 
coordination number z = 4. 
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FIG. 3. The decomposite-spin representation of the 

Heisenberg model (J and infinite Jf) and the next-nearest 
(Ji) and third-neighboring interactions (J2) are shown for 
the pairwise {k a ,kh) tetrahedrons. 



Since the sign of the coupling constant of the XY 
model on a bipartite lattice can be converted by the 
7r-rotation around the z-axis concerning one of the two 
sublattices, in the following we discuss the case with the 
negative interaction. For the spin-1/2 case with N sites, 
it is exactly shown that the ground state is unique with 
Sf 0t = by using the Perron- Frobenius theorem |L(| ■ As 
one can imagine naively, the state with Stot = N/2 and 
Sf ot = is a good zeroth order approximation. In fact, 
the Gutzwiller-type wavefunction constructed from this 
state is shown to be an extremely good variational state 
p| , ^2[ . When applying their results for our z — 4 case, 
the square of the long-range order (LRO) and the short- 
range magnetic correlation are given by; 



(M x2 )/N 2 = 2 15 • 3 4 • 7~ 8 = 0.4604 • 



r k+5, 
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= -1/7, 



(7) 
(8) 



where M = J2k=i ^ k an< ^ ^k+s means the nearest neigh- 
bor pseudospin next to the ffc pseudospin. In the XY 
model, the Hamiltonian has the U(l) symmetry about 
the rotation around the z axis in the spin space. There- 
fore the direction of the order parameter in the t x -t v 
plane cannot be determined even for the infinite system. 
The U(l) symmetric general order parameters, M(6), are 
defined by the ground-state expectation value of the fol- 
lowing operator; 



M (0) = (M x cos 9, M v sin 9, 0) , 



(9) 



where (M x ) = (M v ) — M > 0. Concerning every local 
tetrahedron, M(9) belongs to the E representation of X^. 
Although the effective Hamiltonian, Eq. (Q), has the 
continuous U(l) symmetry as represented by Eq. (Bh, the 
original model possesses only the finite group (T^) sym- 
metry. The U(l) symmetry is fictitious one valid only in 
the lowest-order treatments with the pseudo-spin repre- 
sentation. Suppose that we take account of the higher 
order interactions like biquadratic term and others, then 
this spurious U(l) symmetry should be broken. In order 
to understand how this lowering of the symmetry acts, 
we make use of the Ginzburg-Landau type argument by 
expanding the free energy as a function of the order pa- 
rameter M(9), actually by its magnitude M and phase 9, 
just around the second-order transition point. The U(l) 
symmetry, obtained by the perturbation, corresponds to 
the second-order expansion of the free energy by M(9). 
This is because, on every tetrahedron, the product rep- 
resentation E x E (= A\ + A2 + E) contains the unique 
totally symmetric one (j4i), which is a constant indepen- 
dent of 9. Now we proceed to higher order corrections. 
Since under the time-reversal symmetry only t z changes 
the sign, while both t x and t v remain unchanged, the 
third-order expansion gives the next leading term, intro- 
ducing the anisotropy to fix 9. By calculating the third 



order invariant, which is obtained by reducing the prod- 
uct representation Ex Ex E, we find that the anisotropy 
is proportional to cos 36*. Thus we get a set of the sta- 
ble states, {M (tt), M(±tt/3)} or {M(0), M(±2tt/3)}, de- 
pending on the sign of the M 3 term. These two sets only 
differ in the overall sign and are essentially the same. 
Let us illustrate the symmetry property of the 



order parameter by taking M(0) 



Z^fc=i< 



^2k = i(\uk)(uk\ — \vk)(vk\) as an example. The non- 
vanishing long-range correlation between rj? and t£, re- 
sults in the different population of the \u) and \v) states 
on a tetrahedron. Therefore the parity symmetry con- 
cerning the bonds vertical to the c-axis is broken in the 
ground state, since the \u) and \v) states are character- 
ized by the odd and even parity concerning these bonds, 
respectively. Three different values of 9 in each set of or- 
der parameters correspond to the three equivalent choice 
of the cubic principal axis; M(27r/3) and M(— 27r/3) 
characterizes the break-down of the parity symmetry con- 
cerning the bonds vertical to b and c axis, respectively. 
These three states form a domain structure in a real sys- 
tem. 

Now we extend the present results for positive Jp also 
and summarize the ground-state properties of Eq. (|6|) by 
taking c-axis as a principal axis (6 — 0, n). Depending 
on the sign of J' = J 2 — J\, the uniform (J' < 0) or 
the staggered signed (J' > 0) summation of the (rf ) all 
over the tetrahedrons defines the order parameter. In the 
classical picture, for negative J' the ground state shows 
the ferromagnetic LRO of the transverse components of 
the spin chirality \u) (\v)) when 8 = (tt), while for 
positive J' the antiferromagnetic LRO of the \u) and \v) 
states. For a local tetrahedron problem the compressed 
and elongated Q v mode along the c-axis stabilize the \u) 
and \v) states, respectively 0. Therefore, for the ferro- 
magnetic J' model the order parameter, M(0) (M(ir)), 
induces the uniform Q v lattice distortion with c/a < 1 
(> 1). On the other hand, for the AF J' the staggered Q v 
distortion is induced. Accordingly, under the assumption 
that the inter-tetrahedron interactions are more impor- 
tant than the local magneto-elastic couplings, the inter- 
tetrahedron couplings (Ji and J 2 ) lift the many- fold de- 
generacy of the spin-singlet manifold by the second-order 
phase transition. However, introduction of the magneto- 
elastic coupling, even for a weak case, will induce a small 
first order structural distortion, leading eventually to a 
weak first order phase transition. 

In conclusion, we have studied the ground-state prop- 
erties of the spin-1 pyrochlore antiferromagnet by us- 
ing the tetrahedron-unit decomposition of the pyrochlore 
lattice. In this approach, we assume that the spin- 
singlet manifold spanned by the VBS-type wavefunction 
well-describes the low-energy physics and the essential 
point is how the degeneracy is lifted. We have inves- 



tigated the case where this degeneracy is lifted by the 
inter-tetrahedron interactions caused by the next nearest 
neighbor and the third neighbor interaction, which pro- 
duce the the pairwise effective Hamiltonian of the XY 
type between the pseudospin operators describing the 
tetrahedron singlets. To break the spurious U(l) symme- 
try of the effective model, we have considered the higher 
order anisotropic term based on the symmetry property 
(E) of the order parameter. It has turned out that the 
parity-broken ground state emerges through the second- 
order phase transition. 

Lastly we briefly comment on the relation between the 
present results and the numerical results for the pure 
spin-1 Heisenberg model j|. The singlet ground states of 
a single spin-1 tetrahedron problem are three dimensional 
with Ai+ E irreducible representations. When we define 
these spin-singlet states, with Si + 02 = 03 + S4 = 0, 1, 
and 2, by |a), |6), and |c), respectively, the orthonormal 
basis of the E representation are given by {|C7}, ]V)} = 
{\b), (-2|o) + V5|c))/3}. Considering the fact that the 
symmetrization on every vertex does not affect the lo- 
cal symmetry properties, the parity-broken symmetry of 
the ground state, obtained by using the decomposite-spin 
representation, may correspond to the \U) and \V) states 
in the spin-1 picture. Since Koga et. al suggested the pos- 
sibility of the new spin-gap ground state characterized b 
the J U) state around the isotropic Heisenberg limit 
our result seems to be consistent with their result. 
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